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We explore the possibility of testing modified gravity exhibiting the Vainshtein mechanism against 
observations of cluster lensing. We work in the most general scalar-tensor theory with second-order 
field equations (Horndeski's theory), and derive static and spherically symmetric solutions, for which 
the scalar field is screened below a certain radius. It is found that the essential structure of the 
problem in the most general case can be captured by the program of classifying Vainshtein solutions 
out of different solutions to a quintic equation, as has been performed in the context of massive 
gravity. The key effect on gravitational lensing is that the second derivative of the scalar field 
can substantially be large at the transition from screened to unscreened regions, leaving a dip in 
the convergence. This allows us to put observational constraints on parameters characterizing the 
general scalar-tensor modification of gravity. We demonstrate how this occurs in massive gravity as 
an example, and discuss its observational signatures in cluster lensing. 
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I. INTRODUCTION 



The discovery of the accelerated expansion of the Uni- 
verse [l|, [|| has generated high interest in consistent long- 
distance modification of gravity. Although gravity might 
be modified so as to account for cosmic acceleration, 
modification must be suppressed at short distances be- 
cause general relativity is confirmed to very high ac- 
curacy in the solar-system. Thus, for a modified the- 
ory of gravity to be consistent, it is important to equip 
a mechanism to screen the effect of modification that 
could otherwise persist down to small scales. A wide 
class of modified gravity theories can be (at least effec- 
tively) described in terms of a scalar-tensor theory, and 
the screening mechanisms can be understood by inspect- 
ing the structure of the scalar-tensor Lagrangian. Mainly 
two mechanisms have been proposed so far. One is the 
chameleon mechanism Q , in which a density-dependent 
effective potential is introduced to make the scalar field 
sufficiently massive in a high-density environment. This 
mechanism is utilized for example in f(R) gravity 4|. 
Another mechanism is known as the Vainshtein effect [5[ , 
which operates in models with nonlinear derivative in- 
teractions. In such models, nonlinearity becomes large 
in the vicinity of a matter source to suppress the scalar- 
mediated force. This mechanism is incorporated for ex- 
ample into the DGP model @, 0], the Galileon model Q, 
and massive gravity 0, [13] • 

In this paper, we discuss how the Vainshtein effect 
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works in Horndeski's scalar-tensor theory fllT - tl3l ]. Horn- 
deski's theory is the most general theory composed of a 
single scalar field and a metric that gives rise to second- 
order field equations, and therefore includes a large num- 
ber of concrete models of modified gravity as specific 
cases. We present an algebraic master equation that 
determines a gravitational field around any static and 
spherically symmetric distribution of nonrclativistic mat- 
ter. The similar problem of deriving the Vainshtein so- 
lution in the cosmological background was addressed in 
Ref. (l4j . There, however, complete analysis was hin- 
dered by the complexity of the governing equations. In 
the case where the effect of the cosmic expansion may 
be ignored, the relevant equations are simplified so that 
the problem actually reduces to what has been done in 
the context of the Galileon model 0] and massive grav- 
ity [TH, [lBl ■ It is therefore possible to determine com- 
pletely the conditions under which the Vainshtein solu- 
tions are available. Our solution thus obtained is useful 
in testing modified gravity. 

Among various cosmological observations, we are par- 
ticularly interested in cluster lensing, which can provide 
the convergence over a wide range of radii by combi ning 
strong- and weak-lensing measurements (see Refs. \Yi\ - 
[l9j . and references therein). The precise measurement 
of the structure of galaxy clusters offers us an important 
clue to test modified theories of gravity [2(| Hlj . Working 
in Horndeski's theory, one can examine modified gravity 
endowed with the Vainshtein screening mechanism in a 
generic manner. In this paper, we compute the lens- 
ing convergence of galaxy clusters in Horndeski's theory 
based on the general spherically symmetric solution. As 
a specific example, we demonstrate how nonlinear mas- 
sive gravity proposed recently in Ref. [22| is included in 
our analysis, and explore the possibility of detecting the 
signatures of modified gravity in the context of massive 
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gravity. 

This paper is organized as follows: In the next section, 
we define Horndcski's Lagrangian and explain our proce- 
dure to obtain static and spherically symmetric solutions. 
Then, in Sec. IIII1 we derive the Vainshtein solution and 
present the conditions for the existence of the solution. 
We show the relation between Horndcski's theory and the 
decoupling limit of massive gravity in Sec. IIVI In Sec. |V] 
we study the effect of modification of gravity on gravi- 
tational lensing and explore its implications for cluster 
lensing observations. Section |VI] is devoted to conclu- 
sions. 



II. STATIC AND SPHERICALLY SYMMETRIC 
CONFIGURATIONS IN HORNDESKI'S THEORY 

We consider the most general scalar-tensor theory con- 
structed by Horndeski. The action we study is given by 

S = Sgg [g^u, 4>] + S m [g^ v , ipm] ■ (1) 
The first term is composed of the four Lagrangians, 



S, 



GG 



(2) 



where 

c 5 



K(d>,X), (3) 
-G 3 (0,X)rJ0, (4) 
G 4 (0, X)R + G iX [(D0) 2 - (V^V.0) 2 ] , (5) 



G 5 (cj>, X)G^V^V^ - -G 5X [(D0) 3 
-3D(/»(V M V^) 2 + 2(V M V,0) 3 ], 



(0) 



with X := —g fll/ d fl 4>du(j)/2. The second term is the action 
for matter fields. Here we have written the Lagrangians 
in the equivalent form called the generalized Galileon 
rather than the original one. In this paper, we use nota- 
tion such as G4X ■= dG^/dX and G30 := dG^/dcf). 

We are interested in spherical overdensities on subhori- 
zon scales, so that we will neglect the effect of cosmic ex- 
pansion. In Ref. [lH , the Vainshtein mechanism in Horn- 
dcski's theory has been investigated taking into account 
the background evolution of the scalar field, = 0o(*)j 
and it was shown that Newton's "constant" evolves in 
time through its dependence on 0o(i) and Xo(t) := <Pq/2 
(see also Ref. [23| for a similar argument). However, the 
time variation of Newton's "constant" is strongly con- 
strained from experiments pij . This leads us to assume 
that the scalar field has a negligible time dependence at 
least during the relevant period, 0o ~ const and Xq ~ 0. 

Our background solution is thus taken to be 

ds 2 = 'q^dx^dx 1 ', = 0o = const, X = 0. (7) 

In order for the theory to admit this solution, we require 
that K((f> , 0) = and K^fo, 0) = 0. 



Spherically symmetric perturbations produced by a 
nonrelativistic matter lump on top of the background ([7]) 
can be written as 



ds 2 



[1 + 2<f>{r)}dt 2 + [1 - 2^(r)]8 ij dx l dx\ (8) 



= O + <p(r), 



(9) 



where r is the usual radial coordinate, r 2 = x 2 + y 2 + 
z 2 . The time-time component of the gravitational field 
equations is 



G4 



(r 2 *')' 



G 



i<j>- 



(r 2 ^y 



2r 2 



(G4X — G50) 



2r 2 



+G 5X 



W) 



1. 



6r 2 4 t ' 

while the space-space component reduces to 



(10) 



2G 4 (*' - $') - 2G 4 ^' - (G 4X - G 50 ) 



>\-2 



0, (11) 



where a prime stands for differentiation with respect to 
r. Finally, from the scalar-field equation of motion we 
obtain 



(K x — 2G3<p) 2 2(G3x — 3G4<p X ) 2 — 



2G 



2(2^-$)' 

^2 



+ 4(G 4j y — G50) 



iff (*'-$')]' 



-l I L*4XX — ^LrS^X I ^ r ^5X ~2 



(12) 



Here, all the functions in the coefficients are evaluated 
at = 0o an d X = 0. From now on, we will ignore 
the mass term because we focus only on modified 

gravity endowed with the Vainshtein mechanism. 

One sees that Eqs. (fl"U|) and (|T2"j) can be integrated once 
to give algebraic equations for ip', if?'. In doing so it 
is convenient to use the enclosed mass defined as 



M(r) := 4tt / (-T t *) r 2 dr. 
Jo 



(13) 



The resultant equations coincide with those derived from 
taking the limit = 0o = const and X = in the result 
of Ref. [3. 

Let us introduce six dimcnsionlcss parameters, 
£, 77, fi, v 1 a, and /?, as well as the Planck mass Mp\ and a 
new mass scale A, to rearrange and simplify the expres- 
sions. Those dimensionless quantities are related to the 
coefficients in the above equations as 



M 2 

Gi = G40 = Mpi£, 



(14) 
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and 



both of which arc dimcnsionlcss. 



K x - 2G 3 = r), 

JL 

A 3 1 



-G3X + 3G40j>c — -7-3) 



Gix — Gb<t> — 



G4XX — -G^x 



G 



5A' 



Wc also define 



Mp 



A 3 ' 

X s ' 

3Mpi 
A6~ 



/3. 



1 M(r) 

MpiA 3 87T7- 3 ' 



(15) 
(16) 

(17) 

(18) 

(19) 

(20) 



Now the gravitational field equations (fTOj) and (fTTj) re- 
duce to 



Mpi 

A 3 r 
M P] 

A 3 r 



-£2; + /3a; 3 + A(r), 

£.x + era 2 + /3x 3 + A(r). 



(21) 
(22) 



Substituting Eqs. (|22j) and (|2T|) to the scalar-field equa- 
tion of motion (|12p . we arrive at 



>(>, A) := £A(r) + (| + 3£ 2 ) a; + [/x + 6a£ - 3/M(r)] x 2 + [v + 2a 2 + 4/3£) x 3 - 3/?V = 0. 



(23) 



Solving the algebraic equation P(x, A) = for x, one 
obtains the radial profile of x in terms of A = A(r). 
It is then straightforward to determine the two metric 
potentials by using Eqs. (|2Tj) and (f22|) . Note in passing 
that if £ = then we have a trivial solution x(r) = 
which is not interesting. In the rest of the paper we 
therefore assume that £ 7^ 0. 

Since there is a sufficient number of parameters, at this 
stage the coefficients of the polynomial P(x, A) are free 
in principle. However, it is important to note that the 
structure of P(x, A) in the most general case is still essen- 
tially the same as the corresponding equation in massive 
gravity: it is quintic and the matter source term A ap- 
pears only in the zeroth-order and quadratic terms. This 
structure allows us to proceed following closely the pre- 
vious analysis in massive gravity [l5|, [l6[ . 



III. THE QUINTIC AND CUBIC EQUATIONS 

In this section, we solve the equation P(x, A) = for 
A 1 and A< 1, and single out a solution appropriate 
for our current purpose in each domain. We then derive 
the conditions under which the two solutions are matched 
smoothly in an intermediate region. The procedure here 
is basically the same as that of Ref. [Ty] . The region 
far from the source corresponds to A <C 1, while it is 
assumed that A 1 in the vicinity of the source, where 
the Vainshtein mechanism is expected to operate. It is 
therefore appropriate to define the Vainshtein radius ry 
by 

A(rv) = 1. (24) 
In the outer region (A <gC 1), there is always a decaying 



solution, 



x t 



7? + 6£ 2 



(25) 



which is obtained by neglecting the nonlinear terms in 
P(x, A). We are interested only in this solution, because 
the other solutions, if they exist, do not correspond to an 
asymptotically flat spacetime. 

The stability of the solution in the linear regime can 
be studied by using the truncated action 



S 



d 4 an 



Ml, 



l + 2£- 



Mpi 



R + rjX 



(26) 



Working in this action, one can derive the stability condi- 
tion easily in the same way as in the Brans-Dicke theory: 
in the Einstein frame, the kinetic term for small fluctua- 
tions has the right sign provided that 



77 + 6£ 2 > 0. 



(27) 



We require the condition (|27p for the stability of the so- 
lution x ~ Xf. The same condition can also be derived 
from dP(xi,A)/dx > 0. (See Appendix A for further 
details.) 

Let us turn to identifying the desired inner solution. 
The inner solution is different depending on whether f3 7^ 
or f3 = 0, because the structure of P(x, A) for A ^> 1 
is crucially different. 



A. /3^0 

In the inner region (4> 1), Eq. (p?3")) reduces to 

P(x, A) w £A - 3/3 Ax 2 - 3f3 2 x 5 w 0. (28) 
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The behavior of the solution to this equation depends 
on the sign of £/3. 



(i) gfl < 

In this case, the second and the third terms in Eq. ([28 
balance, so that the solution is 



A 

T 



(29) 



For this inner solution we find ty'/r oc A 2 / 3 and <f>'/r oc 
A 1 / 3 , which is not at all standard gravity. We therefore 
discard this possibility. 

(ii) > 

In this case, we have the same solution as above, x 3 ~ 
—A//3, which is to be discarded for the same reason. We 
have another solution, 



x± 




= const, 



(30) 



for which the first and the second terms in Eq. ([28]) 
balance. For this latter solution we have the correct 
Newtonian behavior: ty'/r ~ $'/r oc A. This is therefore 
the Vainshtein solution we are looking for. 

Having thus identified the desired solutions in the inner 
and outer regions, we consider the matching of the two. 
For simplicity we focus on the case of £ > 0. The solution 
matching in the case of £ < can also be done essentially 
in the same way as follows. 

For £ > (and hence (3 > 0), the outer solution x = 
Xf < can be matched only to the inner solution x = 
X- < smoothly. 1 The smooth matching of the two 
solution is possible if and only if P(x, A) = has a single 
root in (x-,0) for any A > 0. Since P{0,A) = £A > 0, 
if 



P(x-,A) < 0, 



(31) 



the intermediate value theorem guarantees the existence 
of at least one root in (z_,0). Note that P(x-,A) does 
not in fact depend on A. Using the intermediate value 
theorem again in (— oo,a;_) and (0,oo), it can be shown 
that one or three roots exist in general in (x-,0). 
Let x* and A* be the solution to 



dP(x* , A* 
dx 



0, 



d 2 P(x* , A* 
dx 2 



0. 



(32) 



If such x» <S {x-, 0) and A* > do not exist, it is obvious 
that there is only a single root in (x-,0). In this case, 



the smooth matching is possible. If such x* and A % exist, 
P(x, A) = would in general have three roots in 0) 
for some interval of A. Since P(x, A) = for 4 > 1 
[Eq. ([2"8])] has a single root in (x-,0), two of the three 
roots in (x—, 0) disappear as A increases. For parameters 
satisfying P(x*,A*) > 0, the solution corresponding to 
Xf disappears as A increases, implying that the smooth 
matching is impossible in this case. However, if 



P(x*,A,) < 0, 



(33) 



the solution corresponding to Xf remains and hence the 
smooth matching is still possible. Indeed, there is a single 
root for any A > in this case because a simple manip- 
ulation shows that P(x e (A), A) increases with increasing 
A, where x e (A) is the locus of extrema in (x-, 0). 

Summarizing, the outer and inner solutions can be 
matched smoothly provided that Eq. ([3"Tj) and cither of 
the following two conditions are satisfied: (i) G (x_, 0) 
and A* > do not exist satisfying Eq. ([3"2"j) ; (ii) such 

G (x-, 0) and A* > exist, but they satisfy Eq. (|33|) . 



B. P = 

Here again we focus on the case of £ > for simplic- 
ity. The problem reduces to solving the cubic equation, 
P(x, A) = 0, where now 



V 
2 
2a 2 



P(x,A) -> £A+ M + 3£ 2 )x+((i + 6aZ)x 



The inner solution is 



x 3 . 



(34) 



(35) 



(We assume that v + 2a 2 ^ 0.) This solution shows the 
correct Newtonian behavior: ty'/r ~ $'/r oc A. 

The inner and outer solutions can be matched 
smoothly only for v + 2a 2 > 0, which is also required 
from stability of the inner solution. Smooth matching 
also requires that there is a single root in x < for any 
A > 0, or, cquivalently, that P(x, A) has no local ex- 
trema in x < 0. The two local extrema are in x > 
provided that 



/2 + 6a£ < 0, 



(36) 



where we used the condition (|27|) . Otherwise, one must 
require that the discriminant of dP(x, A)/ 'dx is negative, 
i.e., 



{y + 2a 2 ) (tj + 6£ 2 ) > - (ji + 6a£) 2 , (J> + 6a£ > 0.(37) 
o 



1 Since dP(x±, A)/dx = — 6/3x±, the solution x = X— is stable, 
but the other one, x = x+, is not. 



Summarizing, smooth matching is possible if Eq. (|36[) or 
Eq. (J37]) is satisfied. 
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IV. DECOUPLING LIMIT OF MASSIVE 
GRAVITY 

Let us confirm that the conditions for smooth match- 
ing indeed reproduce the previous resuit obtained in the 
context of massive gravity [HI, Eft- To do so, we start 
with finding out the concrete form of K, G3 , G4 , G5 cor- 
responding to the decoupling limit of massive gravity. 
The correspondence can be seen more clearly if we move 
to the covariantized version of the decoupling limit La- 
grangian, i.e., the "proxy theory" proposed in Rcf. 

It turns out that the proxy theory corresponds to 



25] 



K = = G 3 
G 5 



A/2 
(j4 = — — 



+ Mpi(/> 



Mi 



pi 



aX, 



3^X. 



A 6 



(38) 



In massive gravity, the strong coupling scale A is given 
by A = (m 2 A/pi) 1 / 3 , where m is the graviton mass. 

Since the proxy theory contains the Riemann dual ten- 
sor while the Lagrangian of the generalized Galileon not, 
one may wonder how the former is included in the latter. 
Actually, G5 cx X corresponds to the term containing 
the Riemann dual tensor in the proxy theory. The easi- 
est way to verify this is to compare the field equations of 
the two theories. 

From Eq. (f38|) one finds 

7] = n = v = 0, £ = 1, a ^ 0, $ ± 0, (39) 

so that the parameter space collapses to a two- 
dimensional space. The inner solution x_ exists only 
for (3 > and is given by x_ = — 1/^/3(3. Let us define 
C := \fP/a. Then, the condition (|3"Tj) reads 



P(x_ , A) = (l - 3 V3C + 6C 2 ) < 0. 

Solving the equation d x P{x*, A t ) — x*9;j:P(x*, A*) 
which does not in fact depend on A*, one finds 



1 x_ 
This exists if 



1 + 2C 2 - (1 + 4C 2 



iiC 4 ) 1/2 



1/2 



(40) 
= 0, 

(41) 
(42) 



The equation P(x, A) = has three roots in (x_,0) for 
some interval of A if 



P(x*,A*)>0 & 0<C< 




0.6. (43) 



Therefore, smooth matching of the asymptotically flat 
solution and the Vainshtein solution is possible provided 
that 



a<0 or ^>Wi_ 
a ~ V 24 



13 



(44) 
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FIG. 1: The profile of function of the radial coordi- 

nate r. The curves are plotted for (a, 0) — (0.5, 0.3) (dot- 
ted red), (0.8,0.34) (dot-dashed green), and (0.985,0.375) 
(dashed blue), respectively. As a halo density profile we 
adopt the NFW model with M vir = 1.34 x 10 15 M Q /h and 
c vir = 13.8. 



Thus, we have confirmed that the previous result flRHtij 
is reproduced. 2 



GRAVITATIONAL LENSING IN MODIFIED 
GRAVITY 



In this section, we are going to relate our spherically 
symmetric solution to gravitational lensing observations. 
To do so, it is instructive to begin with seeing the typical 
behavior of the Vainshtein solution in massive gravity, 
adopting the Navarro-Frenk- White (NFW) halo density 
profile jgl, for the source p(r) := -T t *. (See Ap- 
pendix B for the detailed description of halo density pro- 
files.) Figures [T]and [2] show the profile of x and its deriva- 
tive, respectively, as a function of the radial coordinate r 
for different values of a and /?. The fiducial parameters 
of the NFW model we use are M vir = 1.34 x 10 15 M Q /h 
and c v ir = 13.8 , which correspond to p s = 7.16 x 10 4 p C r,o 
and r s = 145kpc//i, respectively. The strong coupling 
scale is taken to be A 3 = (100i? ) 2 M P1 = (46.4km)" 3 . 
Then, the Vainshtein radius determined from Eq. is 
ry = 209kpc//i. (As the parameters characterizing the 
profile we choose to use the virial cluster mass M v - n and 



2 Note that our notation is different from those in [l5t Il6| . In 
particular, Oours = — Osbisaetai.- 
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r [kpc/h] 9 [arcmin] 



FIG. 2: The radial derivative of a- as a function of r. Param- 
eters and definitions of curves are the same as in Fig. [T] 



the concentration parameter c V i r rather than p s and r s .) 
One can see that x(r) can have a sharp transition from 
outer to inner solutions, depending on the parameters of 
the theory, which leads to a peak in x'(r). This occurs 
at around the Vainshtcin radius. 

Having seen the typical behavior of the radial profile 
x(r), we now move to investigate how the lensing sig- 
nal is modified in massive gravity. We assume that the 
background evolution of the universe does not deviate 
much from conventional cosmology and use the ACDM 
background with f2 m = 0.3, Oa = 0.7, and h = 0.7. 
The background metric (J7]) is understood to define the 
physical coordinates at the location of the lens object. 

The basic quantity in gravitational lensing is the con- 
vergence, k, which is expressed in terms of the sum of 
the two metric potentials <£> + := (<£> + 1 i)/2 as 



(xs - x)x 



A i $j 



As 



(45) 



with x, xg, and Aj_ being the comoving angular diame- 
ter distance, the comoving distance between the observer 
and the source, and the comoving transverse Laplacian, 
respectively. Using the thin lens approximation, we can 
rewrite the convergence as 



(Xs ~ Xl)xi 
Xs 



dx A$ 



+ • 



(46) 



where xl is the comoving distance between the observer 
and the lens object and A is the comoving three dimen- 
sional Laplacian. Let us now introduce a new spatial 
coordinate as Z — &l(x ~ Xl) 5 whose origin is located 
at the center of the lens object. The projected radius is 



FIG. 3: The lensing convergence n as a function of 9 for 
different values of the parameters of the theory. In these plots, 
the NFW profile is used with M vir = 1.34 x I0 15 M e /h and 
c v ir = 13.8. Parameters and definitions of the curves are the 
same as in Fig. [I] The points with the error bars represent the 
observational data for the high-mass cluster AI689 provided 
by Umetsu et al. [H.&Oill. 



written as rj_ = a^XhO, where 9 is the polar angle from 
the axis connecting the observer and the lens object, and 
ol is the scale factor at the lens object. In terms of these, 
the convergence (|4"rJ)) can be written as 



k(6) 



2(XS ~ Xl)xlQl 
Xs 



A 

dZ— <& + (r), 



(47) 



where r = ^Jr\ + Z 2 . Using Eqs. (J2TJ) and <j22j), we find 
that in Horndeski's theory 



1_ d_ 

r 2 dr 



[r 2 <P + (r)] 



A 3 [(ax 2 + 2f3x 3 +2A)r 



31/ 



Mi 



pi 



2r 2 



(48) 



Figure [3] shows the lensing convergence for the NFW 
profile with different choices of the parameters of the the- 
ory. The points with the error bars indicate the observa- 
tional data for the high-mass cluster A1689 provided by 
Umetsu et al. [T9L [28M 3l| . An interesting feature observed 
in Fig. [3] is that a dip appears at a particular polar angle 
corresponding to the Vainshtein radius. The dip is most 
enhanced for the parameters near the boundary of the re- 
gion in which the smooth matching is possible. Clearly, 
this is caused by the sharp peak in x'(r) at the Vain- 
shtein radius, as seen in Fig. [5] We see from Fig. @]that 
the peak location is certainly determined by the Vain- 
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9 [arcmin] 



FIG. 4: The convergence k as a function of 6 for differ- 
ent strong coupling scales. The curves correspond to A 3 = 
(150# ) 2 Mpi (dotted red), A 3 = (100# ) 2 Mpi (dashed blue), 
A 3 = (50i/ ) 2 Mpi (dot-dashed green), and ACDM (black 
solid), respectively. We take a = 0.985 and /3 = 0.375. 



shtein scale. From Fig. |4] we also find that the depth of 
the dip increases as A decreases. 

In Figs. [5] and [5] we compare different assumptions on 
the halo density profile. Two representative profiles are 
considered here: the generalized NFW (gNFW) jsjlp 
and the Einasto |35l - (38j profiles. We see that a dip ap- 
pears at a characteristic polar angle in the gNFW and the 
Einasto profiles as well. The depth of the dip is enhanced 
for larger 7 S and larger T. 

The appearance of a dip is expected to be a generic 
feature of scalar-tensor theories exhibiting the Vainshtein 
mechanism, because the essential structure of the master 
algebraic equation ([25)1 in general cases are the same as in 
massive gravity. This helps us put constraints on scalar- 
tensor modification of gravity through the observations 
of cluster lensing. 

The decoupling limit of massive gravity constitutes a 
subclass with two free parameters in Horndeski's theory, 
which motivated us to use it for an illustrative purpose. 
We would like to point out here that there are several 
caveats to be aware of when putting observational con- 
straints on massive gravity based on our analysis. First, 
there will be some corrections to the decoupling limit 
because m and 1/Mpi are not exactly zero in reality. 
However, the corrections are small enough in the region 
outside the Schwarzschild radius r g of a lens object and 
inside the Compton length of the graviton Thus, the 
decoupling limit can be used safely at around the Vain- 
shtein radius, which is relevant to our purpose, unless the 
graviton mass is so large that the Schwarzschild radius 




9 [arcmin] 



FIG. 5: The convergence k as a function of 9 for the gNFW 
profile with M vir = 1.34 x 10 15 M Q /h, c vir = 13.8 and 71 = 3. 
The curves correspond to 7 S = 0.5 (dotted red) and 7 S = 1.5 
(dot-dashed green), respectively. For comparison, the con- 
vergence for the NFW profile (75, = 1) (with the same the- 
ory parameters) is shown by the blue dashed line. We take 
a = 0.985 and /3 = 0.375. 




9 [arcmin] 



FIG. 6: The convergence k as a function of 9 for the Einasto 
profile with M vir = 1.34 x 10 15 M & /h and r_ 2 = 145kpc//i. 
The curves correspond to T = 0.1 (dotted red) and F = 0.3 
(dot-dashed green), respectively. We take a = 0.985 and 
P = 0.375. 
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coincides with the Compton length. In addition, there 
will be other corrections since we need in fact to consider 
the decoupling limit in the cosmological background. The 
corrections to the background solution ([7]) are expected 
to be Ag^ = 0{H 2 r 2 )r,^ and A0 = 0(M P iH 2 r 2 ), 
where H is the Hubble expansion rate at the location of 
the lens object. 3 This correction can be neglected relative 
to the perturbations provided that r < (r g H~ 2 ) 1 /' i . This 
implies that our analysis can be applied to massive grav- 
ity at least within the Vainshtein radius ~ (rg/n -2 ) 1 / 3 
if m > H . Finally, and most importantly, no sensible 
cosmological solutions have been found so far in massive 
gravity. See Refs. (4~0-|43| for recent attempts to address 
this issue. 



VI. CONCLUSION 



the appearance of the dip in the convergence is not de- 
pendent on the particular density profile and hence our 
result is robust against the different choices of the profile. 
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In this paper, we have considered static and spheri- 
cally symmetric solutions in Horndeski's theory, which 
is the most general scalar-tensor theory having second- 
order field equations. Under the assumption of weak 
gravitational fields, we have shown that the problem re- 
duces to solving an algebraic equation which is at most 
quintic. Interestingly, the structure of the quintic equa- 
tion is essentially the same as the corresponding equation 
derived in the context of massive gravity. By inspecting 
the algebraic equation, we have presented the conditions 
under which the screened solution is obtained inside a 
certain radius in the most general scalar-tensor theory 

Using our static and spherically symmetric solutions, 
we have explored the possibility of testing modified the- 
ories of gravity endowed with the Vainshtein mechanism, 
focusing on cluster lensing observations. For an illus- 
trative purpose and for simplicity, we have considered a 
specific case of Horndeski's theory corresponding to the 
decoupling limit of massive gravity, and discussed its ob- 
servational consequences. The key effect on gravitational 
lensing is that the second radial derivative of the scalar 
field can be substantially large at the transition radius 
from screened to unscreened regions. This results in a dip 
in the convergence, which will be a marker of the Vain- 
shtein scale. We have found that this dip is enhanced 
near the boundary of the parameter region that allows 
for the screened solution. Such a feature enables us to 
put constraints on modified gravity. 

The following simplifications have been made in this 
paper: the effect of the cosmic expansion is neglected 
and the density profile is well described by the conven- 
tional NFW, gNFW, and Einasto profiles even in mod- 
ified gravity. However, we would like to emphasize that 



3 In the acceleratin g b ranch, the correction to the helicity-0 mode 
is of order A 3 r 2 l39l |4(J| . We do not consider this case here 
because the graviton mass should be small m < Ho so that the 
expansion rate does not exceed the observed one. 



Appendix A: Stability of radial perturbations 

Let us study fluctuations propagating around a spher- 
ically symmetric background. For simplicity, we restrict 
the analysis to radial modes and abandon the matter 
perturbations: tp(r) — > cp(r) + Sip(t,r), $(r) — > 4>(r) + 
6$(t,r), *(r) -> *(r) + 5V(t,r), and ST" = 0. 

Expanding the action to second order in perturbations 
and eliminating the metric potentials by using the grav- 
itational field equations, we obtain the quadratic action 
for Sip: 



Ssip = 47r / dtdr 



3/C )' rax \ 2 2 dP fx '\ 2 



.(Al) 



Here we neglected the mass term. Instead, the 
above action may be deduced from the linear equa- 
tion of motion, which takes the form [(r 3 /C)' /r 2 ]d 2 Sip — 
(l/2)(dP/dx)S(p" + ■■■ = 0. The coefficient JC is defined 
as 

fC(r) = \(v + ^ 2 ) + 4aA + 2(^ + 6a£, + 6/3A)x 

+6 {y + 2a 2 + 4/3£) x 2 + 4 (10a/3 + w) x 3 
+300 V, (A2) 

where we introduced a yet another dimcnsionless quan- 
tity za through G5xx(<f>o, 0) = —3m/ A 9 . 

In order to avoid the gradient instability of the radial 
mode, one must require that 



dP 

dx 



> 0. 



(A3) 



To avoid the ghost instability one must require that 
(r 3 /C)' > 0. This condition involves the new quantity 
w, which do not appear in characterizing the static and 
spherically symmetric solution. Therefore, the proper- 
ties of the solutions discussed in the main text are not 
directly affected by this requirement. 
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Appendix B: Halo profiles 

In this appendix we summarize the halo profile models 
adopted in the main text. We briefly explain the three 
representativeprofiles: NFW [H, HzJ , gNFW [HHH] , 
and Einasto |35l - l38j profiles. To characterize these den- 
sity profiles, we introduce the radius r_ 2 , at which the 
logarithmic slope of the density is — 2 , and the virial 
cluster mass M v - n = 4-7T Jq V ' t p(r)r 2 dr , which can be de- 
scribed by 

47T 

M vir = —rl ir A vivPc ,(z h ), (Bl) 

where r v - lr is the virial radius, A vrr is the virial overden- 
sity, and p c1: (zi l ) is the critical density at , where zl de- 
notes the redshift of the lens object. We take A vrr = 120 
and z L = 0.183 for A1689. 



2. gNFW profile 

The generalization of the NFW model may be written 
in the form jalQ] 

p(r) = . (B3) 

The NFW profile is recovered for 7 S = 1 and 71 = 3. We 
refer to the profile given by Eq. (|B3|) to the generalized 
NFW (gNFW) profile. For the gNFW profile, r_ 2 = 
r s (2 — 7s)/ (71 — 2) and the corresponding concentration 
parameter is given by c_ 2 = ?" vrr /r_ 2 = c vir (7i — 2)/(2 — 
7 S ). We can specify the gNFW profile with the virial 
cluster mass M v { r , the concentration parameter c_ 2 , and 
the slope indices (75 ,71). 



1. NFW profile 



3. Einasto profile 



The Navarro-Frenk- White (NFW) profile is given 
by @,i 



p(r) 



(r/r s )(l + r/r s y' 



(B2) 



where p s is the characteristic density, and r s is the char- 
acteristic radius at which the slope of the density profile 
changes. For the NFW profile, r_ 2 = r s . It is useful to 
introduce the index of degree of concentration, so-called 
concentration parameter, c vlI = r v - lr /r s . The virial clus- 
ter mass M v ; r and the concentration parameter c v i r can 
be used as the parameters of the NFW profile instead of 
p s and r s . 



The Einasto profile is given by [35H3 



P(r) = p- 2 cxp -- 



- 1 



(B4) 



We can specify the Einasto profile with the virial cluster 
mass M v i r , the special radius r_ 2 , the slope index T. 
Based on the Millennium simulation 44j, the authors of 
Ref. (37| claimed that CDM halos can be more properly 



described by the Einasto profile than by the NFW profile. 
They also argued that the best-fit value of r increases 
gradually with the increase of the virial mass, from T ~ 
0.16 for galaxy halos to T ~ 0.3 for the most massive 
clusters. 
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